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We consider a symplectic extrapolation of the Hubbard model of A'^ fold replicated electrons and 
solve this model exactly in two special cases, at TV = oo in the bosonic sector and for any A'^ on a 
dimer of two points. At A'^ = oo we find a multiplet of collective modes that contains neutral spin 
^ . fluctuations and charged pair fluctuations that are degenerate with each other at zero doping. Our 

solution of the symplectic model on a dimer of two points for any interpolates smoothly between 
A'' = 1 and N — oo without any visible discontinuity. These results suggest that the inclusion of 
charged pairing modes in weakly doped antiferromagnets is essential and that an expansion about 
the N = oo limit is appropriate in this context. 
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I. INTRODUCTION AND MOTIVATION 



Calculational schemes that are perturbative in the Coulomb interaction are known to be inadequate for describing 
strongly interacting electrons in systems such as, for example, transition metal oxides, organic conductors and quantum 
dots and this necessitates the use of non perturbative methods, such as the renormalization group, dynamical mean 

■ field theory and various numerical methods. The 1/N expansion for fold replicated degrees of freedom (with 
] subsequent extrapolation to the physical value of = 2 or A'^ = 3) was successfully applied to magnetic impurities 

: 1^ and in other areas of physics Q|. In almost all applications to electronic systems, this expansion uses semiclassical 
saddle point integration and its success in the case of the repulsive Hubbard model ||] is open to question because 
it is unclear where a suitable saddle point exists. An alternative version of the 1/A^ expansion was discovered by 

' 't Hooft in the context of SU{N) Quantum Chromodynamics and relies instead on dominance of diagrams of 
^-H ' simple topology for N ^ oo. Recently it was noticed |^ that 't Hooft 's topological classification of diagrams can 

■ be applied to the Hubbard model, with collective excitations playing the role that gauge bosons play in Quantum 
Chromodynamics. This raises the challenge of reconstructing the large N extrapolation for electron systems from the 
very beginning in a way that avoids saddle point integration. 

Ch ' The symplectic extrapolation of the Hubbard model 0, @ is interesting because, unlike the SU{N) extension, it 
, conserves the notion of two component spins and permits pairing. In the present paper, we solve this model exactly 
(— I in two special cases (i) in the bosonic sector on an infinite lattice at N = oo and (ii) on a dimer of two points for 
Q , any N. Our solution at A^ = oo shows the existence not only of neutral spin fluctuations but also of charged pair 
O ■ fluctuations that are in the same multiplet and exactly degenerate with the spin fluctuations at zero doping. Our 
solution of the dimer for any A^ indicates a smooth interpolation between A^ = 1 and N = oo without any apparent 
singularity. Taken together these results indicate that (i) both spin and pair fluctuations must be included from the 
very beginning in any theory of weakly doped antiferromagnets and (ii) that the expansion in powers of 1/A^ about 
an appropriate N = oo limit is a legitimate procedure in this context. 

The symplectic extrapolation of Hubbard's model that we consider here is based on U{N,q), the group of 
unitary transformations on A^ quaternions and it leads to the following Hamiltonian Q : 

^U{N,q) = ^,y] {{iJxeipx)'^ {ipxSlpx)] +{i't'^xy +^txyi^t^V (1) 

iV t J symm 

X ^ x.y 

1pxS1px= ^ 1pxai{i(Ty)^pi^xf}i, i^ti^x= X! i^tati^xai 

i=l..JV,a/3=l,2 i=1..7V,Q=l,2 
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Here {}sj„„m denotes symmetrization and the electrons ipxai occur in copies with i = 1..N but we will often 
suppress this extra index. For our purposes invariance of Hij(^pf^q-j under SU{N) transformations is sufficient and 
^x£'4'x and in eq(|]) are indeed invariant under SU{N) because ^xu and ^x2i transform, by fiat, according to 

mutually conjugate SU{N) transformations. If we imagine, for the moment, that the interactions in eq([^) come about 
via exchange of suitable bosons, then both neutral and charged bosons must be exchanged. It will be shown below 
that this naive expectation is correct and that the symplectic extrapolation of the Hubbard model contains indeed 
both neutral and charged collective excitations. 

This paper is organized as follows. In section 2 we solve the SU{N) Hubbard model first for two points and any 
iV using SU{2) pseudo spins and we also find the SU{oo) Hamiltonian in the bosonic sector for an infinite lattice. In 
section 3 we determine the N oo asymptotics of the symplectic model after mapping it onto a model of the double 
exchange type and we also solve the symplectic model on a dimer for any N using the group SP{4) ~ iS'0(5) as a 
classifying symmetry. In section 4 we identify the spectrum of the asymptotic SU (oo) and U{oo,q) Hamiltonians. 
Section 5 serves to counterbalance our extensive use of ancient boson operator methods by reformulating our problem 
in standard path integral language. Our conclusions are given in section 6. 



II. EXACT RESULTS FOR LIMITING CASES OF THE SU(N) HUBBARD MODEL 
A. The SU(N) Hubbard model for a dimer 

For simplicity, we begin with an SU (N) extrapolation ||^ of the Hubbard model |jl^ : 

Hsu{N) = ^Yj^^^ + ^^v^ti'v (2) 
which we simplify even further by speciahzing it to a dimer of only two points : 

hI^^] - ^ im' + n,') + \ Y: i^U2. + ^i^H (3) 

1=1. .Af 

Here Ua = X]i=i N'^ai^ai, a = 1,2 and i — 1..N enumerates the electron species and we have chosen hopping 
amplitudes with a (naive) bandwidth of 1. This Hamiltonian describes electrons on a ladder of length N with hopping 
along the rungs i or tunnelling between two mesoscopic quantum dots on a semiconducting substrate, each doped 
with N electron donors. The position a = 1, 2 on the dimer can be viewed as a pseudo spinor index and H^^^^ can 
be expressed entirely in terms of pseudo spin operators as follows: 

Hi^^^^^T! + Tx + const 

^"=1 E ^«^a/3V'/3. (4) 

i=l..N,a = l,2 

where the constant has physical meaning because it depends on electron filling. Because the Hamiltonian decomposes 
into blocks of distinct pseudo spins we may diagonalize it separately in each block and determine the optimal value 
of the pseudo spin at the end (it is impossible to diagonalize the SU (N) dimer for large N by brute force because the 
dimension of its Hilbert space explodes like 4^). It it is useful to diagonalize the hopping term by interchanging the 
conventional representations of and Tx'. 

Hi^N) ^ + const (5) 

Because s-OasiV— >ooa naive first approximation to the ground state is \s, s^ = — s >. To understand the action 
of the Hamiltonian on states in the vicinity of |s, — s >, we express it in terms of Holstein-Primakoff |12| oscillators: 



r_|_ = a^v 2s — a+a, T_ = v 2s — a+aa, — —2s + a^a (6) 
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and then diagonalize it with the help of a Bogoljubov transformation: 



H = \ll + ^b+b + const + 



(7) 



;7s 



b = a cosh 9 + a smh 9, tanh 26 — — ^^^7— 

1 + — 

,„ , tanh^^ I 

\G> = exp( —a+a+)\Q > 

The use of the Holstein-Primakoff expansion for iV — s- cx) is justified, a posteriori, because the number of a^a+ pairs 
converges in the ground state. The optimal value of the total spin s is the maximal one that is permitted at a given 
fiUing e.g. s = f at half filling 

By numerically diagonalizing the SU{N) dimer for finite N the convergence to its = 00 limit can be checked. We 
consider the gap An{U) between the ground state and the first excited state and which, according to eq(0) should 
converge (at half filling) to \/l + U and our calculation confirms this convergence. Figure 1 displays the gap Ajv(t^) 
for N — 2... 5 and shows that the convergence with N is non uniform in U. This reflects the failure of the limits 
N 00 and J7 — > 00 to commute that follows from eq(||). We therefore expect the extrapolation to the physical value 

= 2 to run into difficulties if the Coulomb interaction is too large compared with the bandwidth. 



gap 




_j — , — , — , — . — , — , — , — . — , — , — , — . — , — , — , — . — , — , — , — ^ u 
2 4 6 8 10 

FIG. 1. Oscillating convergence in A'^ of the gaps A]v((7) of the SU(N) dimer for = 2. .5 towards Aoo{U). The highest 
curve is for N — 2, the lowest for = 3 and the asymptotic one is dotted. 



B. Hamiltonian of the SU(oo) Hubbard model on any lattice 



We now generalize the operator representation of the SU{N) dimer to an arbitrary lattice. Because we want to use 
the notion of a Fermi surface, we transform the SU (N) Hamiltonian of eq(^) to plane waves: 



H 



U 



(8) 



p—q—r—s 



with ij^x = "^7=^ S ^^''^V'pj Bpq = ippipq and Cp = J^^xoe^^'^. To find the 1/A^ expansion of the Bpq operators we 
examine the commutators of their fluctuating components: 



[Bpq^Brs] — SqrBps SpsBqr 

— Sqr {Bps — NSpsUp) — 5ps [Bqr — NSqrUq) + NSpsSqr {Up — Uq) 

Up = < tpptpp >= Fermi function 



(9) 
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We change notation and use subtracted bilinears Bpq = ippipq ~ NSpqUp from now on. The preceding equation shows 
that the operators Bpq for p, q on opposite sides of the Fermi surface behave like harmonic oscillators: 



Bp.q ^ VNbp.q , Bp^q ^ VNb+^ (10) 
with [bp^q,b±J = Sp^-Sq^s 

with, so far, unknown corrections. Above we denoted a level by p when it is empty and by p when it is full. Fortunately 
for us, the expansion, to all orders, of bilinears of fermions in terms of bosons was previously worked out in nuclear 
physics, see Q for a review. Although we will not need this representation in its most general form we quote it from 
ref , according to which the leading density operators are represented by 



Bp - q - 



(VF^) bp,r (11) 
V / q.r 



A ^ h- 

^pq — ^-^q^,p 

and where the square root \/N — A is to be interpreted as a series in By contrast, the non leading density 

operators, with momenta on the same side of the Fermi surface, are simple quadratic forms in the harmonic oscillators: 

^pq ~ ^ ^t,p^,q ~ ^qp (1^) 

r 

Bp.-q^N-Y^b + ^bp,, 



The representation of the algebra (|9|) via eqs(|ll|,y_2|) constitutes a generalization of the Holstein-Primakoff repre- 
sentation of eq(^), but the proof that it represents the fermion bilinears correctly is much more difficult for this 
generalization than for the Holstein-Primakoff representation [p^ . Here we only use the lowest order terms of this 
expansion and we retain as dynamical variables the bilinears Bpq and Bpq with p, q on opposite sides of the Fermi 
surface where they have, from eq(^, harmonic oscillator like commutators 

Up + nq = I : [Bpq.Brs] = NdpsSqs [Up - Hq) + 0(1) (13) 

The non leading Bpq for p, q on the same side of the Fermi surface are expressed, using eqs(|l^, ^l]) in terms of the 
leading operators as follows 



(14) 



Bpq = J^'^i^r - rip) BprBrq + ( ) + "^""^^ 



(our notation no longer distinguishes between filled and empty levels) . We may now use eq(^^ to express the kinetic 
energy in eq(^) in terms of leading density operators 

epBpp = ^ Y I^P - e?! {5p9^9p}.j,™,„ + const + O (j^j (15) 

where we have antisymmetrized Cp with respect to the replacement p q. Using eq([T5|) the Hamiltonian in eq(|^) 
reads as follows: 

^^^W = 2A^ ^ ^P«^- + ^ E |ep-e,|{Sp,S,p}„ + consi-|-0(l) (16) 

p-q=r-s n(p)-|-n(q) = l 

As we shall see later in more detail, this describes RPA like particle hole scattering, provided we restrict the summation 
in the interaction term to n{j)) -\- n(q) = 1 and n{r) + n(s) = 1 so that B^qB^s ^ N . 

The operators B^qB^s in eq(p^) contain, however, also non leading contributions with either p, q or r, s on the same 
side of the Fermi surface, or both. By using eqdl^) to express these non leading contributions in terms of the leading 
operators, we find, to next order in the following three point interaction: 
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p— g— r— s;i 

The interaction H^^^ is '-^ N~^/^ because it contains three operators ~ y/N each with arguments on mutually opposite 
sides of the Fermi surface and it gives rise to a correction of the energy. In this paper we will not make use of 
such non leading interactions. 

Finally we check whether the RPA like expression of eq(]l6|) reduces to our previous result for the dimer. We choose 
a naive bandwidth of one so that Ictt — eo| = 1 and take into account vol = 2 to obtain: 

Hli}fS') = J {bo. + btf + 6+ + const + 0(1) (18) 

It is reassuring that this agrees with our previous result (|^) that was obtained directly via the Holstein Primakoff 
representation of pseudo spin operators. 



III. EXACT RESULTS FOR THE SYMPLECTIC HUBBARD MODEL 



A. Map of the symplectic model onto a double exchange model 

The symplectic model of eq(^ is fairly complicated and to make progress we first map it onto a spin model. A hint 
on the mapping we need is provided by the fact that the operators {ipx^'^'x, {'4'x£'4'x)~^ ■, V'^V'a;} form an SU{2) algebra 
under commutation. This fact becomes obvious after a particle hole transformation ip^i V'xi, "0x2 V'J^2 
down spins: 

J V HV'i^^.i - V'i2^.2) / \sij 

We see that this particle hole transformation |l7| replaces the difficult charge ±2 pairing operators by while density 
fluctuations get mapped onto S^. Applying this map on the symplectic Hamiltonian of eq(^ we find: 

Hu(N,q) HdE ^ -^Y~^ X - ^ X + Y txyflJx<^3^y (20) 
X x.y 

While in the symplectic model of cq(^ up and down electrons transform according to mutually conjugate representa- 
tions under SU{N) they transform according to the same SU{N) representation in the model of cq(po|) and therefore 
all four bilinears ^^^a^yP become legitimate SU{N) invariant operators. The Hamiltonian Hoe is of the double 
exchange type, with only a single band and a Hund coupling of the "wrong sign" that favors antiferromagnetic order 
and we will refer to it as "double exchange model" for short. The chemical potential i^ti^x of the symplectic model 
transforms into ip^asipx in the double exchange model where it acts like a magnetic field that reduces the symmetry 
to spin rotations about the z axis H . 



B. Asymptotic Hamiltonian of the double exchange model 

We first rewrite the double exchange model ( pO| ) in terms of plane waves as 

Hde = X! ~^Pi ' '^''^ epPpp (21) 

p+r—q+s p 
'^pq = V'pa^a/sV'g/J, Ppq = ^pa^qa 

where 'r'afs are the Pauli matrices. The fermion bilinears Bpa^qp = tppa'4'qf3 of the double exchange model now turn 
into RPA like collective modes in the same way as the spinless ip ti^o pa irs did in the SU{N) model of the last section, 
except for replacing p ^ ip,ce), g — > (g, f3) and Up — > npa in eqs(p_J 
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rlpa +nqi3 = 1 : [V'paV'g/3, V'^^V'ss] = NSprSgsSasSff-y {Upa " "g/j) + 0{l) (22) 
npa = ?^g/3 : Bpa,ql3 = ^ X! ^'^''^ ~ '^P") Bp^^r-yBr-y + *^ ( ) + '^O'^'^* (^3) 

The Fermi factors in the double exchange model depend on the spin direction for non zero doping. This will break 
its SU{2) * C/(l) symmetry and split the multiplet of modes. To simplify our analysis, we therefore limit ourselves 
here to zero doping for which the Fermi surfaces in the double exchange model loose their spin dependence. We may 
expand the collective modes tppQ.ipqP of 6q(^2|) over the Pauli matrices ~t' and decompose them into scalar ppq and 
vector 'ff'pq operators according to 

^pc.%P = l E '^Pq^'c^ r^^il^r) (24) 

p=0..3 
'^pq ~ {Ppqi~^ pq) 



Using eqs(22,2J) the modes ppq, 'o'pq can be shown to be independent degrees of freedom, to leading order in N: 

H,q. <s] = 'iNSf.JpsSqr {Up - nq) + 0(l) (25) 

Use of this representation on eq(|2^) for Hde gives us the RPA dynamics of pp^q and "ff^p.g in the imdoped DE model: 



Hde = JSJ^yol E ~^P<i ■ "^^^ + ^ E ~ "^9) ('^9 ~ "-P) [PpqPqp + ~^pq ' ~^ qp\ (26) 



+const + O 



p+s=q+t 
1 



N 



where we antisymmetrized Cp with respect to interchange p q. H^e contains also nonlinear interactions in 'o^pq ■ ~^st 
with either p, q or r, s on the same side of the Fermi surface, but we will not write them down here. From our previous 
discussion we know that the operators 5*^ , 5*^ correspond, respectively, to chargedand neutral modes and, similarly, 
a*pq decomposes into charged (cr^^) and neutral (Cp^) bosons. The basic reason for the appearance of charged collective 
modes in the symplectic model of eqdlJ) is the interaction of pairs in its Hamiltonian. 



C. Solution, for arbitrary A^, of the symplectic model on a dimer of two points. 



We consider the dimer both in its symplectic and its double exchange version by specializing eqs(p],po[) to two points: 



^dimer ^ ^ ^ \{{i^,e^,)+ (V^.eV^.)}^ 



(27) 



dimer _ 4[/ 



DE 



For a dimer of two points, the Hamiltonian of eq(G6h reduces to 



H 



dimer 

DE 



U 



1 



^ ( CT^Ott + ~(^Tro) + ^I'S'r " ^ol [POttPttO + O^Ott CT^tto] + COnst + 0(— ) 



(28) 



According to eq(|2^) the operators ctq^ = {/Oottj^Ott} may be represented as CTq^ ~ \/2Nb'^ to leading order and 
therefore the symplectic dimer with naive bandwidth lejr — eo| = 1 at = cx3 is represented in terms of harmonic 
oscillators as follows: 



lim H' 

AT— too 



dimer _ ^ 



DE 



+ {'S' ~t + b+bo 



(29) 



By referring to eqs(H, ^ for the SU{N) dimer we can read off the frequencies of the oscillators after Bogoljubov 
transformation and we find luq — 1, tJL.a = vl + 2[/. 

Straightforward numerical diagonalization of the dimer is impossible because its number of states grows as 16^ 
with A^. However, according to exact diagonalization for A^ = 1,2,3, the low lying states of the dimer are included 
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in a much smaller subspace of SU{N) invariant states. To further simplify the diagonalization within the SU{N) 
invariant subspace, it is useful to recognize that the building blocks of the symplectic dimer ipa£'4'b, V'^V'b for a, 6 = 1, 2 
in eq(|2^) generate the algebra of S'P(4) ~ SO{5). This can be seen by choosing two commuting operators from this 
set such as V'l'V'i ^^^d tpt'^^ ^^'^ determining the Cartan weights under commutation with the remaining operators 
and by plotting the two dimensional weight diagram. Opening a textbook on group theory and comparing with the 
Cartan weight diagrams of the classical groups of low rank one concludes that the algebra is S'P(4) ^ SO{5). 

The irreducible representations of SO{5) have been classified with respect to their SU{2) x SU{2) content and we 
only need the special irreducible representation that contains the vacuum state of the double exchange model and 
which is characterized hy Si — S2 — ^, k = 0...N. The kinetic term in this representation can then be determined 
from the reduced matrix elements given in . 



gap 



1 




J . . ^ U 

12 3 
FIG. 2. Monotonic descent and convergence of the gaps Ajv(f/) of the symplectic dimer for N=1..5 towards Aoo{U). The 
highest curve is for A*' = 1, the lowest curve represents N = 00. 

From numerical diagonalization of the resulting representation of the Hamiltonian that we consider at zero doping 
we find the first excited state of the dimer to be a degenerate vector triplet. In order for U to have its conventional 
meaning at = 1 we replaced U ^ ^ in all formulas and plotted the energy An{U) of the vector multiplet relative 

to the ground state in figure 2 for N — 1..5 together with its asymptotic limit Aoc{U) — y/l + U/3. Comparing 
figures 1 and 2 we conclude that convergence of the gaps towards their asymptotic limit is slower for the symplectic 
extrapolation than it was for SU{N). 



IV. THE SPECTRUM OF THE SU{oo) AND U{oo, Q) HAMILTONIANS 

Although eqs( |l6| , ^6|) give the Hamiltonians of the SU (00) and U{oo, q) models in the sector of even Fermion number 
in terms of elementary harmonic oscillators the spectra of these Hamiltonians are not obvious at first sight. This is 
due to a serious overcounting of degrees of freedom, with roughly one boson for any pair of momenta or any pair of 
points. These essentially non local bosons can describe local physics only by virtue of the interaction ~ U being local. 

We wish to determine the spectra of the asymptotic Hamiltonians via the poles of some correlation functions and 
begin with the correlator < Bp-^q-^B^^^^ > of the SU{oo) theory for simplicity. As a further simplification, we adopt 
the more flexible functional integral representation of the partition function which is fairly obvious here because we 
are dealing with a collection of independent oscillators: 

Bpg = VNbpg (30) 

Up + Hq = I : [bpq, 6+ ] SprSqs {Up - Uq) + 0{l/N) 

(we have rescaled the oscillators to have conventional commutation relations). The Lagrangian we need is just a slight 
generalization of the conventional one: 
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L = l^{np- nq) b+gdtbpg - hsu(N) [b, 6+) (31) 



2 

P,Q 



= J DBe-^ with S = dtL = b+Mb 



Correlators (in imaginary time) of harmonic variables such as b and 6+ are well known to be given by the inverse of 
the quadratic form in the imaginary time action: 

< bp,,, {t)bl^^ (0) >=< p^qi ; t| ^1^292; > (32) 

At this point we invoke the locality of the interaction ~ U which, in x space and in terms of the rescaled variables 
b, reads: 

/•/3 rP 

/ dthint = U / dty^btx^xx = U / dty2KvPxy,x'y'bx'y' (33) 

Jo Jo ^ Jo ^ 

with Pxy,x'y'bx'y' ~ ^xybxx 

The operator P picks out the coincident piece of the wave function and satisfies the projector relation = P. We 
decompose the full action into its "free" and "interacting" part by writing b'^Mb = b'^Lb — Ub'^Pb with 



b+Lb = / dt 



Pi9 n{p)+n{q)=l 



(34) 



where L describes free propagation of particle hole pairs with U = 0. According to general principles, the free 
propagator < 6(t)6+(0) > is just the inverse of L : 

1 Tl — Tl 

< 6p,(t)6+ (0) >i7=o= SprSgs < pq;t\y\rs;0 >= Sprdgs^ — — - (35) 

^ Ot — (Ep — e-q) 

The last expression resembles the particle hole propagator that enters into the calculation of bubbles and Lindhard's 
function and this indicates that we are on the right track. Next we switch on U and exploit = P, the locality of 
the interaction, in a way reminiscent of the treatment of a Friedel resonance: 

< b{t)b+{Q) >u^o= ^ = YZJfp (36) 

= Z + ^1^137^4 

P i_^pip P represents propagation and rescattering of particle hole pairs or a sum of bubbles. The preceding 
argument shows that the spectrum of the SU{oo) Hamiltonian corresponds to zeroes of the operator 

l-UPjp) = S^yd{t) - U < xt\j\yO > (37) 

^ / tx,Oy ^ 

The locality of the interaction ~ U leads to a propagator of coincident pairs and to functions of a single label. Fourier 
transforming we find 

< xt\j\t0 >=< 6^^(t)6+ (0) > (38) 

= ^ np-Tlg ^ _J_ ^ i(p-g)x ic^t - 

y2Z^e dt-(ep-eg) PV^^^ i^-{ep-eg) 

p,q ^ ^ ^' P.f?.^ 

But the correlation < bxx{i)bQQ{Q) > for free fields is also calculable via Wicks theorem: 

< 6xx(t)&i,(0) > +0{^) =< (V'^^o)o >no internal label^ x{x) = -G{x)G{-x) (39) 



with G{x) = -< >= -jr^p = ^ E 



gi(wt+px) 

Q + yiop ^ py - e(p) 
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where h'^°'P is the hopping matrix. Using conventional Matsubara techniques it is easy to check that 



We conclude that the spectrum of the boson Hamiltonian corresponds to the singularities of 

1 1 



5,yS{t) -U < xt\l-\yO > S^yd{t) - Ux{^,t) 



(41) 



where x{^) (V'^V'a;) (V'o^V'o) > is a susceptibility. At half filling, in d—2 dimensions and at Q = (tt, tt) the Fourier 
transform of x('af,i) is well known to diverge logarithmically to +oo because of an incipient antiferromagnetic 
instability. We conclude that the asymptotic Hamiltonian -ffst/(oo) describes bubbles that have a pole near Q — (tt, tt) 
at half filling. This pole in the spin fluctuation channel is also seen in neutron scattering on cuprates p^ . 

At half filling the spectrum of the asymptotic Hamiltonian i?(7(oo,ij) eq(p6|) on an infinite lattice is very similar 
to that of the SU{oo) model. Again the Hamiltonian describes bubbles and again the spectrum is dominated by 
fluctuations at Q = (vr, tt) except that there is now an extra band index. The excitations form a triplet of spin and 
pair fluctuations that are exactly degenerate at half filling. 



V. EXTRAPOLATING FROM N ^ oo TO N = 1 



The aim of the present paper was to establish exact results on the symplectic model in order to see whether the 
expansion in powers of about A^ = oo in this model makes any sense. The boson expansion helped us find the 
asymptotic Hamiltonian for a finite number of points but as we have seen in the last section this method becomes 
inconvenient for dealing with infinite lattices where the bosons seriously overcount the degrees of freedom. In this 
section we return to main stream methods of condensed matter physics and discuss how to extrapolate to from N = oo 
to N =1. 

From the preceding sections it follows that a treatment of doped antiferromagnets and the Hubbard model at low 
temperature near half filling requires a self consistent treatment of electron propagation and spin and pair fluctuations. 
This may be done either using the topological expansion given in ||] or by using the more conventional saddle point 
method. For simplicity we adopt the saddle point method, although the topological expansion is presumably the 
more powerful approach. We may set up a 1/N expansion of the double exchange model by conventional Hubbard- 
Stratonovich transformation: 

[/ /-> \2 1 

D(l)exp-— \^(p ^ - I's'xj = const, ^ -tfj^iaapi^pi (42) 

Z = J D(j)DiPexp-S, S = J dtL 

X 

where t^°^ is the hopping amplitude. Like in preceding work on the SU{N) Hubbard model at large N ( [||) the 
phase factor "i" is essential to get the correct coefficient +U for the coupling in the double exchange model and again 
the factor "i" causes no difficulty because (j) is purely imaginary at the saddle point. To see this, we integrate over 
the fermionic degrees of freedom and find the stationary point of the effective action:: 



Z = I D(j)DTPexp-S ^ I Dcjjexp-Seff (43) 



N 



J J2 ~^ldt - log det (dt + t^°P + fic73 - iVu ■ 



const 



SSeff = ^i(f>\x,t) = IVUTr-a" < xt\ — —\xt >= ^ < > 

(there are A^ species that contribute to < S*" >, hence the factor 1/A^). So the hypothesis that fcj) is real is consistent, 
but difficulties may be anticipated for /i = i^/Uc/y^. Numerical calculations indicate that the spin fluctuations are 
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quickly killed by doping while the pair fluctuation channel remains intact. The likely explanation is that the spin 
fluctuations correspond to longitudinal fluctuations of (j) the 3 direction while the pair fluctuations correspond 
Goldstone like rotations about the z axis. 

One may also expand Sef / to second order in (j) and conflrm the existence of degenerate spin and pair fluctuations 
at half filling and N = oo without ever using boson operators (but we do not see how to avoid them in the dimer 
problem) . A detailed study of the extrapolation to = 1 via coupled integral equations will be given in a separate 
paper. 



VI. CONCLUSIONS 



In this paper we have derived two exact results that serve as first steps towards a 1 /N expansion of the symplectic 
Hubbard model. Firstly, we diagonalized the symplectic dimer for any N using an underlying 5*^(4) ^ SO{5) 
classifying symmetry. We found a smooth interpolation, without any singularities, between = 1 and N — oo which 
suggests that the expansion, in powers of 1/A^, about the point A^ = cx) is indeed legitimate. 

Secondly, we determined the Hamiltonian of the symplectic model for N = oo and found a multiplet of collective 
modes containing charged pair fiuctuations and neutral spin fluctuations that are exactly degenerate with each other 
at half flUing. This tells us that pair fluctuations in doped antiferromagnets are essential degrees of freedom 

On the way to establishing these results we found analogous results for the SU (N) extrapolation of the Hubbard 
model and we also indicated how Holstein-Primakoff like boson operators may be replaced by more familiar and 
more powerful Green's function techniques. The map of the original symplectic model onto a magnetic model that 
is reminiscent of a double exchange model proved crucial to our arguments. This map is interesting in its own right 
because it may provide us with a simple physical picture of the superconducting instability. 

The extrapolation from A^ = oo to A^ = 1 requires a self consistent treatment of fermions interacting with their 
own triplet of collective modes and taking into account mutual renormalization of bosonic and fermionic fluctuations 
in a way similar to the FLEX approach. In the light of our results, the FLEX approach is more closely related to the 
SU{N) extrapolation of the Hubbard model than the U{N,q) one as it ignores the pairing fluctuations of the latter 
model. The integral equations of the extrapolation to A^ = 1 follow essentially from the topological classification given 
previously in but the details remain to be worked out. 

Several scenarios for the physics of pair fluctuations have been suggested, ranging from preformed pairs to 
postulating a quantum phase transition at a critical doping Q. Also, various authors have added pair fluctuations 
on top of antiferromagnetic ones and obtained interesting results It is fair to say, however, that the details of 
the competition between spin and pairing fluctuations still remain to be understood. 
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